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Abstract 

We extend the definition of algebraic entropy to semi-discrete (difference-differential) equa- 
tions. Calculating the entropy for a number of integrable and non integrable systems, we show 
that its vanishing is a characteristic feature of integrability for this type of equations. 

1 Introduction 

> 

Algebraic entropy pQ was introduced as a measure of the complexity of discrete systems, originally 
' for maps over finite dimensional spaces and finite order discrete (difference) equations. It relates 
to general notions of complexity of maps [2Tll3|. Its use as an integrability detector for such 
systems is by now firmly established |14H17| . Its definition has been further extended to infinite- 
dimensional systems, that is to say partial difference equations [T8l[l9]. and it provides a very good 
characterisation of integrability for these systems as well [20J . In all cases integrability causes a drop 
in complexity, making the entropy vanish, while the entropy is nonzero for generic non integrable 
\ systems. 

We give in Section [2] a generalisation of the concept to systems which mix difference and 
differential equations. These semi-discrete equations constantly appear, for example from partial 
discretisation of continuous systems or in the description of symmetries of integrable equations, 
and the "low growth requirement" has already been used for such systems, in conjunction with 
singularity confinement in the earlier Ref. }21] . In Section[3]we recall the link between the singularity 
structure and the value of the entropy. In Section [1] we describe how to calculate the entropy. 

In Section [5] we calculate the entropy for two discretisations of the nonlinear Schrodinger equa- 
tion: the integrable one given by the Ablowitz-Ladik system [22j and a more naive non integrable 
one. In Section [6] we calculate the entropy of Yamilov's form [23] of the discrete Krichever-Novikov 
equation (see also Ref. [S]). In Section [7J we apply the definition to nonautonomous systems, the 
master symmetries of the Volterra, modified Volterra, and discrete Calogero-Degasperis equations 
given in Ref. [25] . 

The outcome of all our calculations is that a vanishing entropy is, here again, a characteristic 
feature of integrability. 

We conclude by indicating some directions for further studies. 
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2 Definition 



We deal with difference-differential equations of the form 

A + Bu n - k {t) 

Un+lit) = c + D Un . ki ty (1) 

where t is a z^-tuple of "times" [t±,t2, ■ ■ ■ , t u ] and A, B, C, D are rational functions of u n (t), . . . , 
u n -k+i(t) and a finite number of their derivatives d tq u n ,d? qtr u n , . . . ,d tq u n _ k+1 ,d? qtr ii n - k+1 , . . . , 
with respect to t\, t2, ■ ■ ■ , t v . The functions A, B, C, D may explicitly depend on the times t\, . . . , t u 
and on the discrete index n (non autonomous case). 

Such equations define invertible recurrences of order k + 1, that is to say a very large number 
of interesting equations. 

The first step is to define from Eq. ([1]) a homogeneous map by projectivisation. The simplest 
presentation of this process is to write 



U n = Xl/Xk+1 

U n -1 = X2/X k +1 

Un-k+1 = Xk/Xk+1 
Equation ([TJ) defines ci map from {Un > ^n — 1 •>••••> ^n— fe+1 

) to (un+i, u n , • • • , Un-k+2) which trans- 
lates into a map </? from the (k + l)-tuple x = [Xi,X2, ■ ■ ■ , X k+ i] to its image y = [Yi, I2, . . . , ife+i]. 
All components of y are polynomials in X\,X2, ■ ■ ■ , X k +i and a finite number of their derivatives 
with respect to the times t%, £2, • ■ • 5 1„. If we attribute to all X's and any of their derivatives a 
weight 1, and weight to all coefficients, even if they depend on the times or the discrete indices, 
these polynomials are homogeneous with some weight w\. 

We then iterate the map <p. Since we have homogeneised the coordinates, any factor common to 
the k + 1 components should be factored out. Doing so we generate a sequence of (k + l)-tuples, all 
with a well defined weight (using the simple fact that derivatives with respect to the times satisfy 
Leibnitz' rule). This defines a sequence of weights Wi,W2, 

The sequence {w m } satisfies the property 

wi +m < wi ■ w m . (2) 
Definition: The entropy is defined as 

e = lim — log(io m ). (3) 

m— >oo 771 

Claim: The previous property Eq. ([2]) of the sequence of weights ensures the existence of this 
limit. 



and 



U n +1 
Un 



Yl/Y k+1 

Y 2 /Y k+1 



Un-k+2 



Y k /Y k 



+1 



3 Singularity structure and entropy 

There is a deep link between the singularity structure [26] of the evolution and the entropy. The 
defining equation ([1]) is such that the map ip has a polynomial inverse ip. This allows us to define 
two polynomials Ku> and K4, as in Ref. [T]: since ip and ip are inverse of each other, their product 
appears as a mere multiplication of the components by some polynomial 

(ip ■ p){x) = K ip (x).id(x) and (ip ■ ip)(x) = K^(x).id(x), 
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where id is the identity map. 

During the iteration of ip, the appearance of some factor F common to all components means 
that the hypersurface given by F = is sent by the corresponding iterate of <p on a singular set of 
(p (singularity means that the polynomial image is [0, 0, . . . , 0]). The only possible factors are thus 
the factors composing and their proper images by ip. We will give some examples of this in the 
sequel. 

If there is no factorisation, the sequence of weights is purely exponential and 



Claim: Any drop in degree during the iteration comes from the possible presence of common 
factors, and is related to the singularity structure for the geometrical reason explained above. 
We will illustrate this fundamental property in some of the examples studied below. 

4 Actual calculations 

One may of course calculate the iterates of tp. The main obstacle is the size of the explicit expres- 
sions, which rapidly exceeds the capacity of hand calculations as well as computer-aided ones. 

The first issue is that we work with an algebra generated by an infinite-dimensional set: the 
X's and their derivatives. One should however notice that, at every finite order, one uses only a 
finite subset of the generators, and manipulates only polynomials. This makes the calculation of a 
few of the first terms of the sequence {w m } possible. 

To go further we will use the method described in Ref. [27J. We calculate the iterates of 92 
starting with a rational initial condition, of the form [X±(ti, £2, • • • , t u ), . . . , Xk+i(ti,t2, ■ ■ ■ , t u )], 
with all the X's polynomials in the times appearing in Eq. JT|). 

One then produces the sequence of degrees {d n } of the iterates. The point is that, for generic 
initial conditions, this sequence retains the feature of the sequence of weights {w m } we are interested 
in: its asymptotic behaviour. In some cases, and with a judicious choice of initial conditions, the 
two sequences {d n } and {w m } may coincide exactly. 

The next step is to extract the value of the entropy, which is an asymptotic quantity, from the 
first few terms of the sequence of weights (resp. degrees). The standard heuristic method is to fit 
a generating function of the sequences with rational fractions. 



The method has already been shown to work remarkably well for maps and lattice equations [17020] . 
and leads to extremely simple rational fractions with integer coefficients. This is related to the 
algebraic properties of the entropy, in particular the still conjectured fact that the entropy always 
is the logarithm of an algebraic integer [l,28j, but we will not dwell on that here. 

Once the first terms of the sequence are fitted with the Taylor coefficients of a rational fraction, 
the method becomes predictive, and is a good check of the validity of the generating function. The 
location of the smallest pole of the denominator (Qd(s) or Q w (s)) of the generating function gives 
the exact value of the entropy. 

Remark 1: When the entropy vanishes, the growth of the sequences is polynomial. A quick 
test is to calculate the successive discrete derivatives of the sequences. Linear growth is equivalent 
to a bounded first derivative, quadratic growth is equivalent to a bounded second derivative, and 
so on. 

Remark 2: The exact value of the entropy may sometimes be extracted from the singularity 
analysis. In one of the cases we will be able to do that (see Section [5.2)1 . 



w n = w{, 



e = log(u>i). 




or 
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Remark 3: It is quite possible to perform the calculation of the iterates without going to the 
homogeneous description. The expressions are then rational fractions rather than polynomials, and 
the degree to retain is the maximum degree of their numerators and denominators. The sequence 
of degrees one produces is usually different from the one given by the homogeneous form. It 
nevertheless yields the same value of the entropy, since the discrepancy between the two sequences 
is at most a multiplicative factor smaller than or equal to the order k + 1. 



5 Two semi-discrete non linear Schrodinger equations 

One source of semi-discrete equations is the incomplete discretisation of partial differential equa- 
tions. The nonlinear Schrodinger equation is a typical example of such systems: 

l 1 9 

iut + u xx + 2 1 -ix | u = 0, 

with u a complex function of a space variable x and a time variable t. If only the spatial coordinate 
x is discretised, we get a semi-discrete equation with one continuous time variable (y = 1). We will 
examine two cases. The first is the integrable one given by Ablowitz and Ladik [22j : 

i d t Uj + (uj-i - 2uj + Uj + i)/h 2 + \uj\ 2 (uj-i + Uj+i) = 0. (5) 

The second one is a slightly more naive one 

i dtUj + (uj-i — 2uj + Uj + \)/h 2 + 2\uj\ 2 Uj = 0. (6) 

Here h is the discretisation step. Both cases are perfectly sensible as discretisations of Eq. ([5]). 
However they have been shown to behave very differently [29J, the reason being that only the form 
Eq. ([5]) retains one of the main features of the original equation: its integrability. We will evaluate 
their respective entropies. 



5.1 Ablowitz-Ladik system 

as a system for the real and imaginary parts qj and rj of Uj, with some 



We rewrite equation 
rescaling: 



d t q n = Qn+i -2q n + g n _i + q n r n (q n+ i + q n _i) 
-d t r n = r n+1 -2r n + r n _i + q n r n (r n+1 + r n _i). 



The homogenisation proceeds by setting 



X 



in. — , - 



Y 



In— l — y i 



u 

Tn-l = y- 



The map ip in the homogeneous coordinates [X, Y, Z, U, V] reads 



<P : 



Y 
Z 

u 



( V 2 (X' + 2 X - Z) - X{W + YZ) \ 
V 2 (-Y' + 2 Y - U) + Y(W - XU) 
X(V 2 + XY) 
Y(V 2 + XY) 
V{V 2 + XY) 



which is a degree 3 birational map over the algebra generated by 

X, X',X", . . . , Y, Y', Y", ...Z, Z', Z", . . . , V, V, V" . . . 



(7) 
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where X' = dtX, etc. 
The inverse is 



Y 
Z 
U 



( 



Z(V 2 + UZ) 
U{V 2 + UZ) 
V 2 {Z' + 2 Z - X) - Z{W + UX) 
V 2 



\ 



-U' + 2 U - Y) + U(W - YZ) 
V{V 2 + UZ) ) 



The direct calculation of the iterates of cp produces 

{«;„} = 1, 3, 9, 19 

which is not long enough to guess the asymptotic behaviour. 

The specialisation of the initial condition with linear functions of t, which provides the simplest 
nontrivial rational initial values for q n , r n , produces the sequence 

{d n } = 1, 3, 9, 19, 33, 51, 73, 99, 129, 163, 201, 243, 289, . . . 

The sequence is (redundantly) fitted by the generating function 



9{s) 



1 + 3s z 
(1 - sf 



showing quadratic growth, vanishing entropy, i.e. integrability. 

Consider degree drop and singularity pattern. The first drop in degree comes with ip 3 , since 
3 x 9 — 19 = 8. We know that this happens when <p hits some singular point. The pattern is 
the following: the surface S + : V 2 + XY = is sent (blown down) to [— X 2 , V 2 , 0, 0, 0] by ip. 
At the next step it goes to [0, 0, X 2 , — V 2 , 0], a singular point which tp blows up to the surface 
X_ : V 2 + UZ = 0. The third iterate sends "smoothly" S + into and the drop of the degree is 
due to the presence of the factor E^_. 

5.2 Naive discretisation of the nonlinear Schrodinger equation 

Consider now equation (j6|) in terms of the same variables q n , r n : 

d t q n = Qn+i -2q n + q n ^i + 2q 2 r n 



-d t r n = r n+ i -2r n + r n _i + 2q n r 2 n . 
With the same notations as in the previous case, the map ip reads 

/ 



(8) 



Y 
Z 

u 



zv 2 
uv 2 

(2Z - X)V 2 - 2Z 2 U + Z'V 2 - ZVV 
(2U - Y)V 2 - 2ZU 2 - U'V 2 + UVV 
V 3 



The sequence of degrees we get is purely exponential 

{d n } = 1,3,9,27,81,243,. . . 

It is straightforward to prove that the sequence of weights is exactly the same; that is, there is no 
degree drop, w n = 3 n . 



5 



The inverse of ip is related to ip by a linear similarity. If A is the permutation 

A : [X,Y,Z,U,V] — > [Z,U,X,Y,V], 

then 

■0 = A • tp ■ A. 

The polynomials k^, and k^, are just 

i^ip = K-ij) = V . 

The "hyper surf ace" V = is sent by ip to one of its fixed points [0, 0, Z, U, 0]. Consequently the 
later iterates cannot hit any singular point, and there will be no factorisation. This proves that the 
entropy is exactly e = log (3) in contrast with the vanishing entropy of the integrable discretisation 
of Ref. [22]. 



6 Generalised symmetries of integrable equations 

A source of integrable semi-discrete equations is the symmetry approach to integrability of discrete 
systems |30H36| . A fundamental equation is Yamilov's form of the discrete Krichever-Novikov 
equation |23j. The equation is 

R(u n+ i,u n ,u n -i) 

O t U n = , (9) 

u n +i — U n -l 

with six free parameters fa, i = 1 ... 6: 

R(u n+ l,U n ,U n -l) = (fa Un + 2fau n + fa) U n+ i Un-i + (fau n 2 + fau n + fa) (u n+1 + Un-l) 

+fau n 2 + 2fau n + fa. (10) 

Using an obvious notation, the homogeneous version of the corresponding map reads 

-Z [faX 2 Y + faX 2 Z + faXYZ + faZ 2 (Y + 2X) + faZ 3 + Y(X'Z - XZ')\ ' 
X [faX 2 Y + faXZ (2 Y + X) + faYZ 2 + faXZ 2 + faZ 3 - Z(X'Z - XZ')] 
Z [faX 2 Y + faZX (2 Y + X) + faYZ 2 + faXZ 2 + faZ 3 - Z(X'Z - XZ')] 

The sequence of degrees obtained with degree 1 initial conditions is: 

{d n } = 1, 4, 10, 20, 34, 52, 74, 100, 130, 164, 202, 244, 290, 340, 394, 452, . . . 

fitted by the generating function 

:i + s )(i + s 2 ) 









H 


Y 


H 




Kz) 





g(s) 



(1-8 



,3 



meaning quadratic growth and vanishing entropy. 
Consider the singularity structure. We have 

with 

A = Z, B = faX 2 Y + faXZ (2 Y + X) + faYZ 2 + faXZ 2 + faZ 3 + Z (-X'Z + XZ') 

C = 2 XZ [(fa fa + fa fa - fa fa) X 2 + (fa fa + fai - fa\ + fa fa) XZ + (fa fa - fa fa + fa fa) Z 2 

+ (fa fa ~ Pi) X 4 + (fa fa - /3 5 2 ) Z 4 + (XZ' - X'Z) 2 . 
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For example, at step ip 2 the product A 3 B factors out causing a degree drop of 6 = 4x4 — 10 = 
3 + 3. 

At the next step the factor is B 3 C <p*(B) causing a degree drop of 20 = 4x 10 — 20 = 3x3+4+7 
(here (p*(B) = is the equation of the proper image by ^ of the surface of equation B = 0, i.e. of 
the pull-back by <p). So the geometrical interpretation of the degree drops is perfectly in order. 

Remark: Changing in an arbitrary way the expression of R completely changes the result: 
for example introducing a new independent parameter /3 7 as the coefficient of u 2 (ti n _i + u n+ i) 
in Eq. (USD leads to a sequence of degrees 1,4,10,24,58,140,338,816,1970,... , fitted by g{s) = 
(1 + s) 2 /(l — 2s — s 2 ), meaning an entropy of log(l + V2). 

7 Non autonomous systems 

Three interesting semi-discrete equations are provided by master symmetries of the Volterra, mod- 
ified Volterra, and discretised Calogero-Degasperis equations, respectively Eqs. (6), (20) and (26) 
ofRef. [25]. 

• The first equation is 

d t u n = u n [(e + n + 2) u n+ i + u n - (e + n - 1) u n _i] . (11) 

It is nonautonomous in the discrete variable n, but the coefficient e is constant. 
The corresponding map reads 

'X'Z - XZ' - X 2 + (e + n - 1) XY S 
(e + n + 2)X 2 
(e + n + 2) XZ 

where X' = dtX, etc... 
One finds the sequence of degrees 

{d n } = 1, 2, 4, 8, 13, 20, 28, 38, 49, 62, 76, 92, 108, 127, 146, 166, 188, 211, 236, 262, 290, 319, 350, . . 
fitted by the generating function 

(1 + 2S 3)(1- S 12 + S 13) 

91 s ) = 5 • 

(l + S )(l-s) 3 

We have quadratic growth of the degree, vanishing entropy, and integrability. 

The second equation is given by 

d t v n = (A - v n ) v n [(e + n + 1) v n+x - (e + n - 1) u„_i] + A v n (v n - A/2). (12) 

Here the parameter A depends on the time and satisfies the condition 

t \ = A 3 /2. (13) 

Equation (|12p is thus non autonomous both in the discrete and the continuous directions. 
For this case, the map ip is 

2 Z(X' Z - XZ') + 2XY(XZ - X)(e + n - 1) + XXZ(XZ - 2Xf 
2X 2 (XZ - X)(e + n + l) 
2XZ(XZ - X)(e + n + l) 









H 
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We get the sequence 



{d n } = 1, 3, 9, 19, 33, 51, 73, 99, 129, . . . 

which has quadratic growth: the second derivative is constant (value 4). 

Notice that if we do not impose the correct time dependence of A given by Eq. f)13[) . sup- 
posing that A is constant, the entropy does not vanish anymore. We get the sequence 
1, 3, 9, 22, 51, 116, 262, 590, 1327, . . . fitted by the generating function 



9{s) 



l + s 2 -s 4 



(1- s)(s 3 - s 2 -2s + l) 



and giving an entropy of approximately log(2. 246979604) [ logarithm of the inverse of the 
smallest root of (s 3 — s 2 — 2 s + l)]. 



The third equation reads 



4 d t v n = (l-v, 



/& 22W e + n e + n-1 \ 2 
(b - a v n ) + a *v n 



(14) 



V n+ l+V n Vn + Vn-l, 

with a = A + \i and b = A — [a, the two coefficients A and fi verifying the same equation (fl~3jh 
We will not write explicitly the formula for the map 99 here. 
We get the sequence 

{d n } = 1,4,13,28,49,... 
which has constant second derivative (value 6). 

If one took a and b constant, forgetting their time dependence, the sequence of degrees 
would be 1, 4, 15, 42, 107, 264, 643, . . . showing exponential growth, and non vanishing entropy 
log(l + y/2). 



8 Conclusion 

All our calculations, without exception, show that the vanishing of the entropy is the hallmark of 
integrability for the extended class of semi-discrete systems, as it was for the purely discrete ones. 
Our study can be complemented in the future in a number of directions: 

• In order to completely validate the calculations made with degrees rather than weights, it is 
sufficient to examine the factors dropping from the successive images under the iteration, and 
check that all factors are in agreement with the geometrical description given in Section [3l 
This will guarantee the absence of spurious factorisation, due to the possible non-genericity 
of the initial conditions. 

• We have not described here equations with more than one time (y > 2), as for example the 
2-dimensional infinite Toda field theory 

d 2 xt u n = exp (u n+1 -2u n + u n -i) 

which can be turned into an equation of the form (fT]) by the change of variable exp (u n ) = v n . 
Our definition applies directly to this type of higher dimensional systems, and the result is 
the same. 
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• Another direction to explore is lattice equations, that is to say systems with more than one 
discrete index, to which the definition can be adapted. 

• Finally, the homogeneous setting we have introduced is more general than the recurrence 
dl]) we started from. The definition we gave can thus be applied to a much wider class of 
semi-discrete equations. 
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